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Abstract 

We observe that successive applications of known results from the theory of positive systems lead to an efficient 
general algorithm for positive realizations of transfer functions. We give two examples to illustrate the algorithm, one 
of which complements an earlier result of [6] . Finally, we improve a lower-bound of [18] to indicate that the algorithm 
is indeed efficient in general. 



Key words: Positive linear systems, discrete time filtering, positive realizations 



1. Introduction 

Given the transfer function 
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of a discrete time-invariant linear SISO system of 
McMillan degree n, we say that a triple AeR" x ", 
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b, c £ R n is an nth order realization of H{z) if it 
satisfies the condition: 

H(z) = c T (zI — A) _1 b. 

It is known that an nth order realization of H(z) al- 
ways exists (see, e.g. [101 Chapter 9]). In this note, 
however, we are interested in the positive realiza- 
tion problem, i.e. finding A, b, c with nonnega- 
tive entries (and possibly of higher dimension M > 
n). The nonnegativity restriction on the entries of 
A, b, c reflects physical constraints in applications. 
Such positive systems appear, for example, in mod- 
eling of bio-systems, chemical reaction systems, and 
socio-economic systems, as described in detail in 
|10|14|15] . A thorough overview of the positive re- 
alization problem and related results has recently 
been given in 4 , while for a direct application in 
filter-design we refer the reader to [5] . 

The existence problem is to decide for a given 
transfer function whether any positive realization 
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A, b, c of any dimension M exists. It is known that 
the constraint of positivity may force the dimension 
M to be strictly larger than n, see Q], [6], [18] for 
different reasons why this phenomenon may occur. 
The minimality problem is to find the lowest pos- 
sible value of M. These problems have been given 
considerable attention over the past decade. The ex- 
istence problem was completely solved in [2] and [9] , 
cf., |16|17|llj . while a few particular cases of the 
minimality problem were settled in |8|13|19|3|23|22j . 

The state of the art of the theory is therefore 
rather two-sided. On one hand, there exists a gen- 
eral and constructive solution [219] to the existence 
problem which, however, is inefficient in the sense 
that it yields very large dimensions, even in trivial 
cases. On the other hand, the minimality problem 
is solved only for particular classes of transfer func- 
tions, and a general solution seems to be out of reach 
for current methods. 

In this note we first observe (Section [2]) that an 
appropriate combination of known results leads to a 
constructive, efficient, general algorithm to solve the 
existence problem in close-to-minimal dimensions. 
We observe that a repeated application of a lemma of 
Hadjicostis [12] leads to the positive decomposition 
problem which, in turn, may be treated by methods 
of [7|20j . In Section [3] we give two illustrative exam- 
ples. In the first we compare the arising dimension 
to that of the earlier general algorithm of [2] . In the 
second we complement the results of [6] by deter- 
mining the minimal value of M for a class of trans- 
fer functions. Finally, in Section [4] we provide a new 
lower-bound on M, improving a result of [18) . This 
latter contribution is independent of earlier results. 

2. The algorithm 

It is known that a necessary condition for the exis- 
tence of positive realizations is that one of the dom- 
inant poles (i.e. the poles with maximal modulus) 
of H(z) be nonnegative real, and there is no loss 
of generality in assuming that it is located at Ao = 
1, see, e.g. [2]. The transfer function H(z) is called 
primitive if Ao is a unique dominant pole. It is also 
known, see [9], that by the method of down-sampling 
the case of non-primitive transfer functions can be 
traced back to primitive ones. Therefore it is cus- 
tomary to assume that H (z) is a primitive transfer 
function with dominant pole at Xq — 1. We shall also 
assume, for technical simplicity, that Ao = 1 is a sim- 
ple pole (this makes the calculations less involved; 



we note that the case of a multiple dominant pole 
can be reduced to the simple pole case as in |171 Step 
4]). Without loss of generality we may assume that 
the residue at Ao = 1 is 1 (see e.g. [2]). 

With these normalizing assumptions, the transfer 
function H(z) takes the form 



H(z) = 



1 
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G(z) 
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(2.1) 



where the poles Xj of G(z) are of modulus strictly 
less than 1, i.e. G(z) is asymptotically stable (note 
that A,s and c^s are possibly complex). 



In the series expansion H(z) = 



k=l 



t k z 



k the co- 



efficients tk are called the impulse response of H(z). 

HH(z) = c T (zI-A)- 1 bthcn t k =c T A fc - 1 bfor 
all k > 0. In particular, tkS must be non-negative 
for H(z) to have a positive realization. We now give 
the main ingredients upon which the algorithm is 
based. The first is the following simple but powerful 
result of Hadjicostis (see |T21 Theorem 5]). 

Lemma 2.1 (Hadjicostis) 

oo 

Let H(z) = ^""^ tjZ ~ J be a rational transfer func- 
j=1 

tion with non- negative impulse response t±, t%, 

For m > 1 let H m {z) denote the transfer function 
corresponding to the shifted sequence t m , t m+ i, . . . , 



Assume that H m (z) 



i.e. H m (z) = y^t m+ j^\z~ 
./ i 

admits a positive realization of some dimension k. 
Then H(z) admits a positive realization of dimen- 
sion k + m — 1 . 

We apply Lemma [2.11 to H(z) as given in (|2.1[) . 
Note that Hi(z) = H(z) by definition, and for each 
m > 2 we have H m (z) = zH m -i{z) ~ t m -\. Hence, 

for eachm > 1, H m (z) = —j+J^ £1=1 {z l%y 
The leading coefficient remains 1, while all other co- 
efficients c^ m — > exponentially as m — > oo (due to 
the asymptotic stability of G(z)). That is, the lead- 
ing coefficient becomes large compared to other co- 
efficients, and this is exactly the familiar situation 
of the positive decomposition problem, which we now 
turn to. 

The task in the positive decomposition problem 
is to decompose an arbitrary transfer function G(z) 
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as the difference G{z) = T^z) - T 2 (z), with Ti(z) 
and T 2 {z) both admitting positive realizations (see 
[7113119120] ). By rescaling, one may assume that 
G(z) is asymptotically stable, and then the usual 
approach is to take a one-dimensional positive sys- 
tem T 2 (z) — where < w < 1 is larger than 
the modulus of any pole of G(z), and R is a suffi- 
ciently large positive number. Then T\(z) — G{z) + 
T 2 (z) can be shown to admit a positive realization 
which, in some cases, turns out to be also minimal 
7; 13 19 20 . For our purposes, the essence of these 
results can be summarized as follows: for any prim- 
itive transfer function, as long as the partial fraction 
coefficient of the dominant pole is significantly larger 
than all other coefficients (as in H m (z) and Ti(z) 
above) there exist efficient methods to construct posi- 
tive realizations. We shall not list all relevant results 
of [7113119120] concerning the positive decomposi- 
tion problem; instead, we give as an example The- 
orem 8 of |7j, which handles all transfer functions 
with simple poles. 

Theorem 2.2 (Benvenuti, Farina & Anderson) 
1 1 "-i 

+0W _ +£^, 

j=l 3 

where G(z) is a strictly proper asymptotically sta- 
ble rational transfer function of order n, with sim- 
ple poles. Let Vj (j > 3) denote the interior of the 
regular polygon in the complex plane with j edges 
located at the j-th roots of unity. Vj can formally 
be defined in polar coordinates as in [J: 

V r .= {{p,6) : p cos(^±^-tf)<cos^, 

forfc = 0,l,...,j-l}. 

Let Ni be the number of non-negative real poles 
with positive residue in G(z) and let N 2 denote the 
number of other real poles in G(z). Let Ns denote 
the number of pairs of complex conjugate poles of 
G(z) belonging to the region V3, and let Nj (j > 
4) denote the number of pairs of complex conjugate 
poles of G(z) belonging to the region 

3-1 

Vj \\Jv m . 

m— 3 

If all CjS are sufficiently small then H{z) admits a 
positive realization of dimension N = (n—l)+N 2 + 

J2ti - 2)Nj = 3Nj. 

3>3 j>l 




■■■ V 3 ■ ■ ■ and V 4 

Figure 1. The sets Vj. 

Remark 2.1 

The dimension N appearing in the theorem is not 
necessarily minimal but it is a good a priori upper 
bound on the order of the realization. Further, by 
carefully analysing the proof in jTjj, the condition on 
the residues may be given explicitely to be ^ \cj\ < 
2~ 5 / 2 ( where the sum runs only over the j 's for which 
Xj is counted in Nk, k > 2). In such case it can 
be computed that H{z) admits a positive realization 
of dimension N = jNj . For the reader's conve- 
J>1 

nience, we have indicated in the appendix 

The proof in J^j is constructive and it gives a so- 
called cone-generated positive realization (see Sec- 
tion 4)- 

This theorem was later improved and generalized 
in various forms, [13l Corollary 2], [19l Corollary 2], 
[20l Theorem 1, Theorem 2]. These papers also pro- 
vide a number of examples where minimality of the 
arising dimension N can be claimed. Finally, a syn- 
thesis of all these results, [201 Theorem4], covers the 
case of H(z) = + G(z) for any asymptotically 
stable rational transfer function G(z) . 

Theorem 2.3 (Matolcsi, Nagy & Szilvdsi) 
IfG(z) is any asymptotically stable rational transfer 
function with poles X\, . . . A r of order mi, . . . m r , 
and if all the partial fraction coefficients of G 
are sufficiently small, then the function H(z) — 
— 3j + G(z) admits a positive realization, the dimen- 
sion of which is given explicitely as a function of 
Ai, . . . A r , mi, . . . m r . 

We remark that in the case of simple poles the 
value of N in Theorem 12.21 is better than in Theo- 
rem [231 All the results above are constructive (see 
[7113119120] ). 

Theorems 12 . 2 1 and 1 2 . 31 yield that the following is a 
general algorithm which terminates in a finite num- 
ber of steps for any given transfer function H(z). 

ALGORITHM: 
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Assume H(z) is given as in 
m := 1 

WHILE t m > DO 

IF assumptions of Theorem = TRUE 
THEN APPLY Theorem O 

APPLY Lemma ED 
ELSE 

IF assumptions of Theorem H2]= TRUE 
THEN APPLY Theorem EH 

APPLY Lemma O 
ELSE 77i := m + 1 
ELSE there is no positive realization of H(z) 

Remark 2.2 

At each step, we may apply to H m (z) other known 
constructions from the literature different than The- 
orems \2~M or \2.3\ We included only these two theo- 
rems to keep the algorithm transparent and because 
they guarantee that the algorithm terminates in a 
finite number of steps. Other important partial re- 
sults which can be incorporated into the algorithm 
are given in i8\3\22\j . 

Remark 2.3 

We acknowledge that this algorithm is merely an 
observation that some earlier results in the litera- 
ture can be combined together. Nevertheless, we find 
it an important observation as it provides a com- 
pletely general algorithm. Previously such a general 
algorithm has only been given in JE[. There are good 
heuristic arguments to believe that the algorithm 
above is efficient in terms of producing small dimen- 
sions, and better than the existing algorithm of f^. 
First, the partial fraction coefficients decay exponen- 
tially, so that only a few iterations are needed before 
Theorem \2.2\ or \2.3\ become applicable, and these the- 
orems already provide minimal or close-to-minimal 
dimensions (see Section [3| for a numerical exam- 
ples). Second, the method of involves the time 
development of an n-dimensional "cube" around the 
vector (1, 1, . . . , 1) and, as such, can only produce 
dimensions larger than 2™ (usually significantly 
larger than that). This fact, however, by no means 
diminishes the theoretical significance of the results 
of J^j which provided the first general solution to the 
existence problem for primitive transfer functions. 



3. Examples 

In this section we give two examples. In the first 
we compare the arising dimension of realization with 
that of the algorithm of [2] . In the second we com- 
plement a result of [6] and determine the minimal 
dimension of positive realizations for a class of trans- 
fer functions. 

First, we note that in the case where there are only 
simple poles, the number of iterations needed may 
be evaluated as follows: write H(z) — Y^j=o z -\ 
with Aq = 1 and |Aj| < 1. A simple computation 



shows that H m (z) = Y^j=o 
if m > 



n max cv 



log max | Xi | 



It follows that, 



< 2- 5 / 2 



then J2 

and Theorem l2 . 2l applies (cf. Remark [2.ip . Moreover, 
it is enough to consider those poles that are not non- 
negative with non- negative residues. 

Example 1. Let 

^) = -L_ + ^) = -L_ + 

0.3331328522 z 2 + 0.1984152016 z + 0.1253986950 
z 3 - 0.69055619 z 2 + 0.80189061 z - 0.38920832 ' 
where t(z) is a low-pass digital Chebyshev filter of 
order 3. The partial fraction decomposition is 



H(z) = 



1 



with 



Ai = 0.07522998673 - 0.8455579204 i, A 2 = Ai 
d =-0.01050864690+0.1411896961?;, c 2 = cT 
A 3 = 0.5400962165 c 3 = 0.3541501460. 

Let us now apply the algorithm. The above compu- 
tation shows that after m — 5 iterations, we may 
apply Theorem l2.2l to H§. As N2 = and as Ai , A2 € 
■p 4 \ V3, A 4 = 1, Nj — for j ^ 4 we obtain a pos- 
itive realization of H 5 (z) of dimension N = 5 (cf. 
Remark l2.ip . A repeated application of Lemma I^TTI 
gives a positive realization of H(z) of dimension 9. 
A more careful examination of the proof of Theorem 
12. 21 would show that H3 already satisfies the require- 
ments so that a 7-dimensional positive realization 
exists. We spare the reader the numerical values. 

While we cannot claim that the arising dimension 
is minimal let us compare it with that of the algo- 
rithm described previously in [2] . In short, a specific 
minimal realization of H(z) is considered and the 



4 



time-evolution of a small 4-dimensional cube around 
the vector [1, 1, 1, 1] T needs to be checked to pro- 
vide a system-invariant polyhedral cone, which then 
leads to a positive realization of H(z), The num- 
ber of edges of the cone equals the dimension of the 
positive realization. If one carries out this procedure 
word-by-word for H(z) above the arising dimension 
is 48. 

Example 2. Consider the family of transfer func- 
tions 

4-(5/2) Ar - 2 3 -5 N ~ 2 



H N (z) = 



1 



(3.1) 



2-1 2-0.4 2-0.2 

as in [6] Example 4] . It is proved in [6] that for any 
N > 4 the minimal dimension of positive realiza- 
tions of H N {z) is at least N. Here we prove that 
an A-dimensional minimal positive realization of 
H N (z) does indeed exist for every N > 4. 

For N — 4 the following 4-dimensional positive 
realization of H A {z) is given in [6, Example 3] and 
it is shown to be minimal. 



and 



A = 



1 



/o 

63 + 4^26 



1 



85 







6 51 
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o 22 + 4 ^o 

85 



(3.2) 



Consider now H N (z), N > 4. It is not difficult 
to see that with the notation of Lemma [2~T1 we have 

ffl(z) = ff 2 (2) = ff^M, • • ■ , ff m (z) = 

ff JV+1 - m (z). Let us stop at m = JV - 3, i.e. at 
H m (z) = H A (z) and make use of the realization 
(|3.2[) of H A (z). Then, the application of Lemma l2.1l 
produces a positive realization of H (z) of order 
4 + (m— 1) = N. Note that for producing this mini- 
mal positive realization we use (|3.2|) instead of The- 
orem H2] (we have, in fact, followed the suggestion 
of Remark |2T2"|) . 

We have concluded that the minimal dimension 
of positive realization of H N (z) is N. Let us now see 
what the word- by- word application of the algorithm 
of Section [3] gives. 

The algorithm terminates when Theorem 12. 21 be- 
comes applicable, i.e. when H m (z) — H°(z), that is 
m = N + 1. Then a 3-dimensional positive realiza- 
tion of H°(z) is constructed, and the application of 



Lemnra l2"Tl produces a positive realization of H N (z) 
of order 3 + (m — 1) = N + 3. 

This example is reassuring in that the application 
of the algorithm of Section [3] produces positive real- 
izations of close-to- minimal order. 

4. Improved lower-bounds 

We saw in Section [3] that the minimal order of 
positive realizations of H N (z) is N. Also, it is easy 
to calculate (see [B]) that the impulse response se- 
quence of H N (z) contains zeros, namely ijv-i = 
tjy = 0. A general lower-bound presented in [18] 
gives that in such case the order M of any positive 



realization satisfies 



A/(M+1) 



- 1 + M 2 > N, i.e. M is 



at least w y =gr . In view of the actual minimal value 
M = N a lower-bound of the order of magnitude N 
is welcome (instead of the order of magnitude y/~N). 
In this section we present such an improvement (but 
we note that while the lower-bound of [TH] is valid 
in general, our improvement is restricted to transfer 
functions with positive real poles, as is the case of 
the example of the previous section) . 

Throughout this section we assume that H(z) is 
a given primitive transfer function of McMillan de- 
gree n with positive real poles, and there exists a pos- 
itive integer ko, such that for the impulse response 
sequence of H(z) we have tk = and tk > for all 
k > kg. This means that H (z) is of the form 



i r n, (i) 

— + TT 3 

— 1 ... , [z — XjY 

1 l—l v J ' 



(4.1) 



where Cj € K, < Xj < 1, and Y^j=i n j =n — 1. 

Let the triple (h, F,g) denote an arbitrary mini- 
mal (n-dimensional) realization of H{z) (for canon- 
ical minimal realizations see e.g. pill). Assume that 
there exists a matrix P of size n x M such that for 
some triple (c, A, b) with nonnegative entries: 



FP = PA, Pb = g, c 1 = h 1 P. 



(4.2) 



There is a well-known geometrical interpretation of 
these equalities. Namely, the columns of matrix P 
represent the edges of a finitely generated cone V in 
R n , such that V is F- invariant, and V lies between 
the reachability cone and the observability cone cor- 
responding to the triple (h,F, g). It is known that 
the triple (c, A, b) provides a positive realization of 
H(z). 
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Definition 4.1 

A triple (c, A, b) which arises in such a manner is 
called a cone-generated realization of H(z). 

It is a basic result in the theory of positive realiza- 
tions that a transfer function H (z) admits positive 
realizations if and only if it admits cone-generated 
realizations (see [H]). Here we present a lower bound 
on the order of cone-generated realizations of H(z). 
For this we shall need the following auxiliary result. 

Lemma 4.1 

Let f : E -> K be defined by 

r 
3=1 

where Ai > A2 > • • • > A r > and p/j\ denotes 
a polynomial (with real coefficients) of degree rij. 

r 

Then f has at most R = Rf := 2_X n o + 1) ~ 1 
pairwise distinct real roots. 



Lemma 14. II is proved by induction on R. 
We are now ready to give an improvement of the 
lower-bound of [THj . 

Theorem 4.2 

Assume that H(z) is a transfer function of McMillan 
degree n, with positive real poles, given as in i|4.f)) 
above. Assume also that there exists a positive inte- 
ger ko, such that for the impulse response sequence 
of H(z) we have t ko — and tk > for all k > ko- 
Then the dimension M of any cone-generated posi- 
tive realization of H(z) satisfies M > -^r. 

Proof. Let the triple (h, F, g) denote a minimal 
(n-dimensional) realization of H(z). Consider any 
cone-generated positive realization (c, A, b) of H(z) 
arising from a matrix P of size n x M, as in !4.2l Let 

denote an arbitrary column of the matrix P, and 
consider the sequence g^' := h T F fc ^ 1 e i > 0. Let 
Pi := [pij] be the nonnegative matrix of size M x 00 
defined by Pij := g^ for 1 < i < M and 1 < j. 

Let K := [fcij] denote the infinite Hankel matrix 
composed of the impulse response sequence of H(z), 
i.e. kij := £j-|_j_i. By assumptions imposed on P 
there exists a matrix Q = [qij] of size ooxM, with 
nonnegative entries, such that QP X = K. This is 
true because the fcth row of K is given by kkj = 
t k+ j-i = h T F J '- 1 (F fe - 1 g), and the vector F* -1 g 
lies inside the cone V by assumption. Thus, it may 
be decomposed as a linear combination of the edges 



ei of V with nonnegative coefficients, and [one choice 
of] these coefficients form the fcth row of the matrix 

Q 

Since (h, F, g) is a minimal realization, for an ar- 
bitrary column of the matrix P, the transfer func- 
tion corresponding to the impulse response sequence 
— h T F fc_1 ei is of the form 



H ie '\z) = 



Ci 



1 



(Note that some coefficients Cj and dS I may be 0.) 
The column e^ of P is called dominant if Cj 7^ 
111 H^(z). Delete the non-dominant rows from the 
matrix Pi and the corresponding columns from the 
matrix Q. The remaining matrices (of sizes Mi x 00 
and 00 x Mi for some Mi < M) are denoted by 
P^ om) and Q( d °™). We see that Q(rf°™)p( dom ) < K 
entrywise. Recall that tk = 0, by assumption on 
the impulse response of H(z). This implies that for 
some dominant index i (1 < i < Mi), gjp = 0. 
Otherwise, ki t k — ffe would be strictly positive in 
the first row of K. Considering the second row of K 

we see that &2,fc -i = ^feo = 0> hence g k( '_ 1 = for 
some dominant index i. By the same argument, for 
every 1 < j < fc , gp = for some dominant index 
i. In other words, each of the first fco columns of the 
matrix p( rfom ) contains a zero, and hence there are 
at least fco zero entries in p( dom \ 

On the other hand, g k ^ is the impulse response of 

00 _ 



H^H^ = §=1 + E-=x TTsU ^y- Thus, 9 t 

Ci + Y^=i P{j)(k)*j> where py) are polynomials of 
degree not exceeding nj — 1, for 1 < j < r. Hence, 
Lemma [47T1 implies that there are at most R = (1 + 
Y]j—i Tij) — 1 = n — 1 zeros in each row of p( dom )_ 
This means that the number of zeros in the matrix 



l(dom) 



is at most Mi(n — 1). Therefore, 



fco < # of zeros in p[ dom > < Mi(n-1) < M(n-l), 



and. hence, M > 



n-i ■ 



□ 



Remark 4.1 

If we apply Theorem \4-S\ to functions H N (z) of Sec- 
tion^ we obtain M > N/2. This is still far from 
the actual minimal value N . However, if there are 
more than 3 poles present in H(z) then the geomet- 
ric arguments of JfJjj are difficult to generalize, while 
Theorem\4-.2\ still applies. 
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Remark 4.2 

As mentioned in the "Open Problems and New Di- 
rections" section of IS it is desirable to have tight 
upper and lower bounds on the minimal order of a 
positive realization in general. Note, however, that 
the results of [61 81 and Theorem \4-.S\ above are all 
based on the assumption that the impulse response 
sequence of H(z) contains at least one 0. The only 
other lower-bound known to us is that of \12}/ which, 
however, does not give any non-trivial estimates for 
transfer functions with nonnegative poles. 

What can be said if the impulse response does not 
contain zeros? Unfortunately, we do not have a gen- 
eral approach to this case. As a first step in this di- 
rection we examined the modified family H N,e {z) — 



4-(5/2) 



+ 



3-5 J 



. „ • ,,o f or small values of e. 

„ -1 2—0.4 2 — 0.2 J J 

Note that the impulse response sequence no longer 
contains zeros. Since the dimension of the system 
is 3, we can use elementary (but tedious) geometric 
arguments to conclude that for small enough e the 
minimal order M of positive realizations of H N ' £ (z) 
still satisfies M > N/2. It is not clear, however, how 
to generalize these arguments to transfer functions 
of higher degree (as in Theorem \4.Hty where the geo- 
metric intuition is missing. Therefore, finding tight 
lower-bounds in the general case remains an open 
problem. 

5. Conclusion 



We have observed that recent results in positive 
system theory can be put together to produce an ef- 
ficient, general algorithm to the positive realization 
problem of transfer functions. We have given two 
examples to illustrate the algorithm. In the first we 
compared the arising dimension of realization with 
that of an earlier general algorithm of (2]. In the 
second we examined a family of transfer functions 
given in [6], and determined the minimal order of 
positive realizations. With respect to the minimal- 
ity problem we have proved a new lower-bound on 
the order of positive realizations of transfer func- 
tions with positive real poles, improving an earlier 
general result of [18] . 

Appendix A. Precisions on the paper [7\ 

In this section, we will show how to obtain a quan- 
titative bound on R in Theorem 1 2 . 21 above . This re- 
sult is not really new as it can be obtained directly 
from the proof of that theorem in [ 7i and some simple 



observations. This appendix is thus only included 

here for the reader's convenience. 

Proposition A.l FA Proposition 7] Let H(z) = 

tjw and assume that pe G V m ■ Then 



r/e 



z—pe 1 - 6 



z—pe 

2/7 



for R > '— there exists A + € R™ xm , b +1 c + e 

~ cos^ + + + + 



such that 



Hx{z) = H(z) 



R 



z - 1 



= cl(zI-A + )-% 



Proof. Let us first consider the Jordan realization of 
H 1 : Ht(z) = c T (zI- A) with 



A : 



fpcosO —psin9 0^ 
p sin 9 p cos 9 
1 





^2?y(cos?9 + sin??)^ 




(l\ 


b = 


2fy(cos — shi'tf) 


, c = 


1 




I R ) 




w 



Next, define K by 

/ /2tt\ /2tt 
1 cos — cos — 2 
\m J \m 

sin ( — ) sin (—2 

\ m J \m 

1 1 1 

and, for a > define 



/2tt 
cos — (m 
V m 

sm — (to — 1) 
\ TO 

1 



D(a) = 



(a o\ 

a 

V 001 / 

It was proved in [7] that the positive cone K, a gen- 
erated by D(a)K is A-invariant. To prove the the- 
orem, it is then enough to prove that we can choose 
a in such a way that 



K a C O := {x e M 3 
and that if R > -^L 



:c T A k x>0, fc = 0,l,...} 

_ then b e K, a . 
For the first one, note that if 

x = A(a cos ip, a sin ip, 1) T S M 3 

then 



7 



c T A fc x 



( p k cos kO -p k $mk6 (A 



= A(1,1,1) 



p k smkO p cos k9 
1 



a cos ip 
a sin if 
1 



= A [ap k (cos(v3 + k6) + sin(</? + + l] . 

It follows that, for a = 1/2, A > 0, p < 1, we always 
have c T >l fc x > 0, in particular, /Ci/ 2 C C 

Finaly, note that the cone K a contains the cone 
{A(r, a) : r < cos A > 0}, in particular S /C Q if 



2 3/2 



?7 7T 

i?a to 



(A.l) 



By taking a = 1/2 we can see that this is the case 



2 5 / 2 T7 

as soon as R > # . 



□ 



If the filter is given and that one seeks for more 
precise estimates, one can slightly improve the result 
by taking a — l/2p and still a bit further if = rix 
for some r (necessarily rational) for which one may 
have cos(jff + kff) + sin(j^ + kff) < k < 1 for all 
integers k > and < j < to. In this case one may 
take a = t^-. an d then check for the smallest R for 

which (cos $ + sin ?9, cos ^ — sin i?) 6 "P,„ . 

In the opposit direction, note that R = 2 7 / 2 r] 
works for all to > 3. 

Let us now show how the estimate on R in Theo- 
rem l2.2l results from this. We will take the notations 
from [7j. Let q be the smallest integer so that all 
complex poles of H are in V q and let us write 

Ni (1) q Ni / 



.7 = 1 Z ' 



A 



=2 i=i 



1 



where cf ] > and A^ > 0, ffj 2) 
sponds to the other real poles and, for i > 3 



(0 



(')p-Wj(i) 



where /?) e 
Now, each 



2; 



Z- p 



e ^\U fc < ffc- 

has a one dimensional one- 



dimensional positive realization A~- , 6^, c~- . 



(2) 



> 



„( 2 ), 



According to [7, Proposition 5], if R 

(2) R {2) 

H- ' + ^f-j- has a two-dimensional positive realiza- 



tion Af\bf\ C] 



(2) 



(the estimate is clear from the 
end of the proof of that proposition in [7] and the 



restriction A < 1 is irrelevant). According to the 
improvement of Proposition 7] given above, if 



Rf > 4?7 



(2) 



then H 



W j_ fl i 



2-1 



(0 f,(i) „(*) 



has a i-dimensional 



positive realization Aj ' ,bj ,Cj , A positive realiza- 
tion dimension N (defined in Theorem l2.2p of H is 
then given by 



1.00 



i°l '■■■'°AtJ ) 



„(!) J 2 ) 



(2) 



' c JVi ' c l ) ■ ' ■ ' C N 2 ' 



' 1 ' ' ■ ■ ' C AT„ 



„(«)\ T 



and A + is given in block-diagonal notation by 
A + = diag(A^ \ . . . , A^ ,A^ \ . . . , , . 



A (i) A (g)\ T 



(i) 



A condition in the theorem is thus that R :— Rj 
satisfies e.g. R > |cf | + 2 7 /^ E? ^ 3 E ^ „W. 

This bound may be improved slightly with the above 
remark. 

Note also that this realization has dimension N% + 
2N 2 + 3N 3 --- + qN q . 
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